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Abstract 



We present technical details of the evaluation of the elastic vrvr scattering am- 
plitude to two loops in chiral perturbation theory. In particular, we elaborate 
on the renormalization procedure at the two-loop order and on the evalua- 
tion of the relevant Feynman diagrams that can all be expressed in terms of 
elementary functions. For the sake of clarity, we discuss these matters both 
in the A-component (j)^ theory (in its symmetric phase) and in chiral pertur- 
bation theory. Estimates for the relevant low-energy constants of 0{p^) are 
presented. Threshold parameters and phase shifts are then calculated for two 
sets of O(p^) coupling constants and compared with experiment. We comment 
on the extraction of threshold parameters from phase shift data. 
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1 Introduction 

Elastic TTTT scattering is a fundamental process for QCD at low energies. It provides 
an ideal testing ground for the mechanism of spontaneous chiral symmetry breaking. 
Since only the (pseudo-)Goldstone bosons of chiral SU{2) are involved, one expects 
the low-energy expansion of the scattering amplitude to converge rather rapidly. 
The systematic procedure for this low-energy expansion is called chiral perturbation 
theory (CHPT) i 0, |, 

In the framework of CHPT, the vrvr scattering amplitude was evaluated at lowest 
order in the chiral expansion in |Q and to next-to-leading order in [0, ^. The 
theoretical developments up to 1994 are summarized in |0, §]. The forthcoming 
experimental improvements concerning vrvr scattering are also discussed there. The 
amplitude to next-to-next~to~leading order was recently calculated in a dispersive 
approach in and in the standard CHPT framework in [jlO|. Other developments 
since 1994 can be found in Refs. |]ll]|-[pO[ . 



Both the explicit calculation of the Feynman diagrams and the renormalization 
procedure needed to evaluate the scattering amplitude at the two-loop order are 
quite involved. The main aim of the present article is to detail the methods used in 
|]T0[| for that purpose. 

The amplitude contains a nontrivial analytical part that can be expressed in 
terms of logarithms that generate the cuts required by unitarity. In addition, there 
is a contribution which consists of a polynomial in the external momenta. Our 
calculation reveals the dependence of the six coefficients in this polynomial on the 
pion mass and on the low-energy constants of both O(p^) and O(p^). We find that 
the complete two-loop amplitude can be expressed in terms of elementary functions. 

We found it useful to first illustrate in Sect. 2 the corresponding calculation 
in the A^-component (f)'^ theory (in its symmetric phase), where the algebraic ma- 
nipulations needed are simple. The loop expansion is presented together with the 
renormalization at both the one- and two-loop levels. A technical issue that also 
comes up in the CHPT calculation is the use of the equations of motion (EOM) 
for the counterterms. We show how different forms of the counterterms at the one- 
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loop level lead to the same generating functional after a proper redefinition of the 
counterterms at the two-loop order. 

In Sect. 3, the analogous calculation is performed in CHPT. The renormaliza- 
tion procedure is again discussed in detail, both for minimal subtraction and for 
the modified minimal subtraction that we actually use. For the effective chiral la- 
grangian, the role of EOM is similar as in (j)^, albeit algebraically more involved. In 
Sect. 4, the results for the pion mass, the pion decay constant and the tttt scattering 
amplitude are presented to 0{p^) in the modified minimal subtraction scheme. Here, 
we also compare our results with the calculation of Ref . ^ . 

In Sect. 5, we perform a numerical analysis to compare with experimental infor- 
mation. For this purpose, we first derive estimates for the low-energy constants of 
0{p^). We include the effects of vector and scalar meson resonances and of K and rj 
contributions. We then discuss numerically the sensitivity of the amplitude to those 
constants. We emphasize that the main uncertainties are due to the couplings of 
O(p^). We use two sets of these couplings to illustrate the uncertainty: The first set 
used in [|lOl is essentially based on phenomenology to O(p^) For the second 



set, we use the present calculation to extract the couplings li, I2 from the observed 
D-wave scattering lengths. For the two sets, we then analyse threshold parameters 
and phase shifts and compare them with available experimental results. 

Sect. 6 contains our conclusions. Some technical aspects on the calculation of 
two-loop diagrams, the off-shell scattering amplitude of and analytic expres- 

sions for the scattering lengths as well as for the coefficients 6j in the 0{p^) amplitude 
are relegated to four appendices. 



2 A^— component (f) to two loops 

As we have mentioned in the introduction, we illustrate the loop expansion in CHPT 
first with the A^-component 0^ theory in the unbroken phase. For this purpose, we 
proceed in a manner as analogous as possible to the chiral expansion. In particular, 
we rely on the external field technique and use a procedure that is scale independent 
at each step, as is done in CHPT. We also comment on the role of counterterms 
that vanish upon use of the EOM and give the relation to the MS scheme. 



2.1 The loop expansion 

The loop expansion is equivalent to an expansion of the generating functional in 
powers of ^ - we therefore explicitly display h in this section. The lagrangian is 

1 n °° 

C = Ud,fd^ct^-M^fcl^)-l{fct^f-ff-Y.n^C^, (2.1) 

where the symbol collects the fields 0i, . . . , ^^r, 

(f)^ = (01, . . . ,0Ar) . 

Below, we will use an analogous notation for any multicomponent vector, e.g., 
= (01,02,03), etc. The quantities denote the counterterms that remove the 
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ultraviolet singularities at order . They are linear combinations of 0(A^)-invariant 
polynomials of dimension < 4. Considering the external field / to be a spurion 0{N) 
vector and using partial integration to eliminate df^(f)'^d'^(f), one is left with 

Pi = \f^<i> , P2 = \m^^<P , P3 = \{<f<P? , Pa = 0^/ . (2.2) 

As is shown below, one may eliminate one of these polynomials by use of the EOM. 
Here we discard P4 and write 

Ci = ga^P , C2 = g'b^P , ... , 

= {Pi,P2,P3). (2.3) 

In the following we use dimensional regularization. As a result, the vectors a and 
b become M-independent functions of the space-time dimension d, divergent as 
(i — ^ 4. The external field / allows one to construct the generating functional by 
evaluating the vacuum-to-vacuum transition amplitude 



exp{iZ[f]/h} = N / [d^] exp{iS/h} , S = 




The normalization constant ensures that Z[0] = 0. Expanding the right-hand 
side in powers of h generates the series 

z = Zo + hZi + r?Z2 + 0{r?) . 

To arrive at the explicit expressions for the components Zj, one considers fluctuations 
around the classical solution that is determined by the external field through the 
EOM, 

(M^ + n)^" + (70^00" + /"^ = ; a = l,...,Ar. (2.5) 
Below, barred quantities always denote quantities evaluated at ( p.5|) , e.g., 

5 = 1 dxL^mj) , (2.6) 

etc. We write the fluctuations in the form 

= + e , (2.7) 
use translation invariance of the measure, [d^] = [rf^], and find 
exp{?Z//i} = N ^^^{iSITi} X 

/ \di\ exp{^ I dx ^-eDi + F(e, m ■ (2.8) 

The differential operator D acts in coordinate and flavour space, 

= (M^ + U)5''^ , a"^ = gW4)5''^ + 2^y^} . (2.9) 
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Furthermore, considering the fluctuations to be of order h^^"^ , the quantity F starts 
at order fc'^'^ . For the evaluation of the generating functional at two-loop order, it 
is sufficient to keep terms of order fi^ in F, 



d^ 



'ab 



i (ai □ + a^M^) + (2.10) 
We then arrive at the following expressions for the components Zj, 

Zo = -\j dx {P4-2P3} , (2.11) 

Zi = dx[i{Hl + D^^a))-2Ci] , (2.12) 

Z2 = -9'ldxdy^l{2G% + G,y{G%)}^y (2.13 a) 

+ ydx{ [(G..)^ + 2{Gl)] - [. . .]^^,} (2.13 b) 

- Y y dxdy(j)l {{G,^,^) + 2G^^ + idi^} G^^y x 

{{Gyy) + 2Gyy + idly} (2. 13 C) 

+ig j dx [{d2yGy^)y=^ - {■ ■ ■)y=^,^=o} (2-13 d) 

-jdxC2. (2.13 e) 
Here we have introduced the propagator in the presence of the external field /, 

i"00 

G^y = / rfA(x|e-"^||/) , 

DfG% = 6'''6^''\x-y) . (2.14) 

The symbol (A) denotes the trace of the matrix A in fiavour space, and (p^ = 0(x). 
The various contributions to Z2 are illustrated in Fig. |l]. 

The numbering is such that the contribution ( |2.13 a|) corresponds to Fig. |l|a, 
etc. We refer in the following to Fig. [l|a (b) as the sunset (butterfly) diagram, 
respectively. Below, we will also need the free propagator, 

= Gxy\f=o , (2-15) 
which is a diagonal N x N matrix in flavour space. 



2.2 Renormalization 

The above expressions have a well-deflned limit as the space-time dimension d 
reaches its physical value d = 4 only for an appropriately chosen singular behaviour 
of the vectors a and 6 as d — > 4. In this subsection, we display and comment these 
singularities and start the discussion with the one-loop graphs. 
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Figure 1: The two-loop graphs in (j)^ theory. Graph a corresponds to the contribu- 
tion (|2.13 a|), and analogously for the diagrams b-e. The solid-dashed lines stand 



for the propagator (|2.14| ). The black circles denote a vertex generated by the la- 
grangian J-'l^,^^ in (p.l|), and the filled squares stand for the contributions from the 
counterterm Ci. Finally, the contributions from the counterterm C2 are indicated 
with an open square in graph e. 



2.2.1 One— loop diagrams 

The singularities at one-loop order may easily be identified by writing 
(ln(l + Do V)) = (Do V) - ^(Do-VDq-V) + ■ ■ ■ , 

where the terms omitted are finite at ci = 4. Evaluating the trace and keeping only 
the singular terms at c? = 4 gives 

2 

1 1 dx (ln(l + Do V)) = -g{N + 2)Tm J dx 0^0 - |-(iV + 8)Tm J dx {$^^f + ■ ■ ■ , 

where the tadpole integral iTm is the free propagator at coinciding arguments in 
coordinate space, 

1 /"oo M'^~'^ 
Tm = - rfA(x|e-^^«|x) = jy^T{l - d/2) , (2.16) 
I Jo (47rj^/^ 



and 



Tm ^ Om^Tm = -^^r(2 - d/2) . (2.17) 
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To determine Ci, we work out the singular part in the integrals Tm and Tm- We 
introduce the scale fj,, 



M'"" = /i"™(M//i)""' ,w = d/2-2 , 



and obtain 



with 



M ( 1 

= — ^/i'" - + a(M//i) - 1 + 6(M//x)w; + 0(w;2) 

= ^{- + a{M/^i) + c{M/^i)w + 0{w'')] , (2.18) 

a(M//i) = ln(MV/i2) _ p/^i) _ 1^(4^) _ (2.19) 



The functions h{M/ fi) and c{M/ fi) will not be needed in explicit form in the follow- 
ing. Since we have introduced the scale /i in such a manner that Tm and Tm are 
scale-independent, we may use scale-independent counterterms as well, 



, 




/i^- 1 


i-iN + 2) 


167r2 


w 


1 


f-(iV + 8) 


167r2 





as = \ ^ +4(/i,^)| , (2.20) 

where the renormalized couplings al{fi,w) are finite at w = 0, with 

L|- + 2t/7^a^(/x,w;) = 2(iV + 2) , (2.21) 



and analogously for 03. With the choice (|2.20 , 2.21|) , the one-loop functional Zi is 
scale independent and finite at d = 4. 

2.2.2 Two-loop diagrams 

In the next step, we determine the counterterms that render the two-loop contri- 
butions finite at = 4. For this purpose, we identify the singular part of the full 
propagator 



Sxy ^xy ~)~ ^x^i 

^xy = dM^^xy ■ 



The remainder R 



G xy Sxy ~\~ Rxy 5 (2.22) 
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is finite at x = y. The decomposition ( p.22| ) generates in the two-loop functional 



Z2 terms of the form R^jR^SjRS^ and S'^. A closer examination reveals that all 
nonlocal singularities proportional to R and to R^ cancel out after putting in the 
values of the from Eq. ( ^.201 ). The divergent contribution to Z2 is therefore 
obtained by replacing everywhere in Z2 the propagator G by its singular part, G 
S. We now sketch how the singular pi eces ar e evaluated, and consider for illustration 
the part proportional to G^y in Eq. ( p.l3 a]) that contributes to the sunset diagram 

Fig. After G ^ 5*, we expand the field (py around y = x and are then left with 
divergent vacuum integrals of the form 

dx {A%; A%A,y) (1; (x - yY^x - yY) . (2.23) 

The integrals proportional to A^^.^^ correspond to the sunset integrals H and H^^ 

discussed in appendix The integrals that contain A.j.y can be evaluated in a 
similar manner. Proceeding in an analogous fashion for the remaining contributions 
to Z2, we finally find that the following counterterms render the two-loop functional 
finite at d = 4: 

,2w ^ 2 



, /.^- V[ (5 + iV)(2 + iV) , (2 + iV) ^^ 

' 16^ j 1 ^^^2 + — ~ — (3-a2(/i,w)-a3(/i,w)) + 62(/i,w) 



l£) I ^^^^ + ^ (22 + 5iV - (8 + N)al{ii, w)) + ^^(/i, w) 



(2.24) 



Remark: In Eq. ( p.24|) , the full contribution a2 3(/U, w) to the singular part in 62,3 
occurs. We could have expanded 3 (/i,w) around w = in ( p.20| ) - then, only 
a2 3(/i,0) would occur in (|2.24| ). In order to keep 02,3 scale independent, additional 
terms of order 0{w) must then be added in ( |2.2CI| ). From ( p.24| ) it is clear that these 
would contribute a local part to Z2 that could then be removed with a redefinition 
of the renormalized couplings bl{fi,w). The definition ( |2.2CI| ) allows one to do all 
this in one go. 



2.2.3 The role of the equations of motion 

We come back to the choice of the counterterms Cj[|. In the generating functional, 
they are evaluated at the solution to the EOM. This amounts to a linear relation 
among the Pi, 

E = 2Pi + 2P2 + 4P3 + A , 

E = . (2.25) 

^ The following discussion is adapted to CHPT, where the higher-order lagrangians £4,Ce, . . . 
take the role of the counterterms Ci, C2, . . . . 
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As we already mentioned, one may therefore remove one of the Pi from the list ( p.2| ). 
We may e.g., ehminate P3 instead of P4 as above, 

C,^C[ = ga'^P' , P'^ = (Pi, P2, P4) . (2.26) 

Using the identity 

^Tp, ^ Jp + a'r^E , 

J = (a;-24,a'2-24,-44) , (2.27) 

it is seen that a different choice of the counterterm at one-loop order amounts to 
adding to the lagrangian a term that vanishes at the solution to the EOM. As we 
now demonstrate, the use of C[ requires that the counterterm at two-loop order 
must be adapted accordingly, 

C2 C'2 = g%'^P . (2.28) 

In C2, we use again P - the freedom in the choice of polynomials occurs at each 
order in the perturbative expansion. We denote by Z[, Z'2 the generating functional 
obtained with ( p.26| , p.28|) . For the one-loop functional to be finite at c? = 4, the 

integrals / dxC^ and / dxCi must have the same singularities. According to the 
relation ( |2.27|) , this requirement is satisfied with 

« = (~y'^2 - y,-^) . (2.29) 

The replacement C\ C[ then amounts to 

Ci^Ci + KgE , (2.30) 

with K = —03/4. This transformation obviously leaves the one-loop functional 
unchanged. Repeating the calculation of Z2 with ( p.30| ), one finds 

Z'2 = Z2 + Jdx (C'2 - C2) + Q , 

where Q is a local action, 

Q = —g'^K J dxe^P , 

e"^ = (2ai + K,2a2 + K,4a3 + 6k) . (2.31) 

In addition to the singularities in Z2, the counterterm has to cancel the ones in 
Q as well. Because the latter is a local action, it can be removed from Z2 altogether 
with 

b' = b + Ke . 

In summary, a different choice of the counterterm at one-loop order amounts to 
adding a term to the lagrangian that vanishes at the solution to the EOM, see Eq. 
( p.30| ). That transformation leaves Zi untouched and changes Z2 by a local action 



Q which can be completely removed by an appropriate choice of the counterterm at 
two-loop order. 




Ci + KgE 
C2 + g'e^^P 



Z[ = Z^,Z'2 = Z2 . (2.32) 
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2.2.4 Connection with the MS scheme 



The above renormalization scheme generates scale-independent quantities at each 
step. In particular, individual Feynman diagrams are independent of fi. The 
Green functions so evaluated at one-loop accuracy contain the finite parameters 
M^, g, ttiifi, 0), 0), . . .0. An analogous procedure is used in chiral perturbation 
theory [g]. We find it useful to compare at this stage this scheme with the MS 
scheme in conventional perturbative calculations. There, one starts from the bare 
lagrangian 

^ = lid,^ld'^<Pn - Mlfys) - ^-fifys? (2.33) 
and sets at one-loop order 

gB = ^i-''"9r{^^)^l-^^^9r{^A + 0{gl)^ , 

Ml = M^M jl - + 0{gl) \ . (2.34) 

The two schemes can easily be related by evaluating, e.g., the physical mass Mp and 
the elastic scattering amplitude to one-loop accuracy. This is done in th e follo wing 
subsection in the present scheme. Repeating that calculation with (p.33| , |2.34|) and 
identifying the one-loop expressions gives 

9r{^A=9{l + ^,+0{g')) , 

M'M = M\l + ^^+0{g^)). 

Here and in the following we use the notation 

< = a[(/^,0) , 6[ = &I(/i,0) . 

From the scale dependence of the renormalized couplings 03 , 03 one reads off the 
scale dependence of M^, g^ that guarantees scale independence of physical quantities. 
Analogous relations hold to all orders in the perturbative expansion. 



2.3 Two— and four— point functions 

For illustration, we evaluate the two-and four-point functions. It allows us to ex- 



plain the methods used in Ref. |]T0| for the evaluation of the tttt scattering amplitude 
in CHPT. 



^Of course, since the present theory is renormaUzable, these parameters occur in such a com- 
bination that they can be lumped into, e.g., the physical mass Mp and the physical coupling 
constant. 
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e ^ 

Figure 2: Diagrams that contribute to the two-point function at two-loop order. 
A filled square (open square) stands for the contributions from the counterterm 

Cl(C2). 




2.3.1 Two— point function 

To evaluate the two-point function 



i{Q\Tr{x)<p\ym,onn = 5"'A' 



xy 



= 5'^'{N-'{A,,y)+0ig)) , (2.35) 

we may either calculate the relevant Feynman diagrams displayed in Fig. ^, or 
determine the contribution of order to the generating functional, because 



(2.36) 



We first proceed in this manner and find 



dX(t)iK^y<i)y + 0{?l\<i)'') , 

K,y = 5\x-y){Any + B) + CA,y + DAly 
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with 

A = l-g\ , 

B = M^-gm^-g^{h2M^ + {N + 2){a^TM + m'tM)} , 
C = g'm' , D = -2g\N + 2) , 

= a2M^ + {N + 2)Tm • (2.37) 
Next we consider the Fourier transform 

M^-p^-B + M^A-2g^{N + 2)H{p^) ' ' ^' ' 

where H{p^) denotes the sunset integral (|A.9| ). Expanding H{p^) around = 
gives at d = 4 

A'(P^) = j^^, + R,{P') + 0{g') , (2.39) 



with 

IGvr 



1 + 9^777^ {{N + 2)(^ + a(M//.)) - 6^) + 0{g') , (2.40) 



and 

R^ip') = 2g'iN + 2)Tip'). (2.41) 

An integral representation for the twice subtracted sunset integral H{p'^) - which 
is regular at p^ = Mp - is given in (|A.12| ). Finally, the physical mass Mp can be 
easily evaluated from Eqs. (|2.38| , |A.13|) to two loops. The expansion starts with 

= + gm^ + 0{g^) 

= M'[l + ^[(AT + 2)(a(M//i) - 1) + a^] + 0(/)} . (2.42) 

We add a remark concerning the calculation of the two-point function from the 
Feynman diagrams displayed in Fig. |^. Instead of evaluating the diagrams di, . . . jdg 
individually, one may calculate the tree and tadpole graphs shown in Fig. where 
the dashed and double lines stand for the propagators 

,J 2 ' ^1 = ^'[1 + 9(^2 + 9\b2 - h)] (dashed line) 
^ M2 = M^ + gm'^ (double-line) 



M| - p2 

and where the vertex denotes the coupling g^l+asg). The sum of these two diagrams 
is 

^2^_^2 ' Mi = M^i+g{l + a3g){N + 2)TM, , (2.43) 

where the argument in the tadpole integral is M2. The expression ( ^.43[ ) is equal to 
the sum of the graphs di, . . . jdg up to terms beyond two loops. 
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+ 



X 



Figure 3: Summing up the graphs di ■ ■ ■ dg in Fig. See text after Eq. (|2.42|) . 



2.3.2 Four— point function 

The number of diagrams at two-loop order becomes quite large - of course even more 
so in the case of CHPT. In the case of the two-point function, we have just seen that 
almost all graphs in Fig. |^ can be summed up by evaluating the tadpole diagram 
shown in Fig. ^. Here we wish to illustrate an analogous method to evaluate the 
elastic scattering amplitude (four-point function) to two loops. Whereas the method 
does not really pay off in the case of the A^-component 0^ theory considered here, 
it turns out to be very useful in CHPT. Furthermore, a similar procedure works in 
the evaluation of form factors^. The four-point function is of the form 

= rr rA#2'^ o-. Tim;ik{s,t,U]pl,pl,pl,pl) , (2.44) 

llil^'^P Pi) 

with 

Pl+P2=P3+P4 , 

S= ipi+ P2f , t = (Pl - P'if , U = {pi- Pif . 

The scattering amplitude is obtained by putting all momenta on the mass shell, 

Tim-ik{s,t,U;pl,pl,pl,pl)\p2=M% = SikSiraA{s,t,u) + CJCI . 

To determine the amplitude A[s,t,u), it suffices to calculate the matrix element for 
the indices 

i = k = l , m = l = N . (2.45) 

The relevant graphs up to 0{g^) are displayed in Fig. ^. The numbers attached 
on the external lines denote the relevant momenta, i ^ Pi. The group indices 
and crossed diagrams are not shown, and mass and counterterm insertions on the 
external lines are not displayed. 

We use the obvious notation 

A{s, t, u) = (7^(1) + (72^(2) + (^3^(3) + 0{g^) . (2.46) 

■^This method has already been used in ^ in connection with the evaluation of the process 
77 — > TTTT to two loops. 
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Ci 



C2 




C3 





C4 



C5 








Figure 4: The elastic 00 00 scattering amplitude to two loops. The numbers 
on the external lines denote the momenta. Crossed diagrams and insertions on 
external lines are not displayed. A filled box (open box) denotes contributions from 
the counterterm Ci(C2). 



The tree-level and one-loop results are 



with 



= -r7r^K + iN + 8)a{M/f^)] 

iOTT 

+ 2{N + 4)J{s) + 4\j{t) + J{u) 



J{z) 



a In h 2 



■ (7 = (1 - AM'^/z) 



(2.47) 



(2.48) 



The calculation at two-loop order may be made more economical by using the 
renormalized one-loop off-shell amplitude as a single (nonlocal) vertex. Then, a 
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Figure 5: Summing up the diagrams Ci ■ ■ ■ C5 in Fig. |[ See text after Eq. (|2.47|) . 



large part of the two-loop diagrams can be obtained automatically by making a 
one-loop calculation where one of the vertices is the one-loop amplitude, and the 
other is the standard tree-level vertex. To be more specific, consider the diagrams 
Fig. Hci,. . . ,C5. They can be calculated by evaluating the integrals indicated in 
Fig. g. The notation is as follows: The double lines denote off-shell legs (the 
corresponding off-shell amplitude is defined in ( |2.44| )). These legs carry momenta I 
and {pi + P2 — I), and group index n. The integration has to be done with weight 

1 



(M2-Z2)(M2-(pi+p2-/)^) ■ 

From this sum, one has to subtract the diagram Fig. ^3. The result equals the 
diagrams Fig. |ci,. . . ,C5. 

One may easily check the correctness of this statement, if one singles out of 
the full one-loop amplitude one particular diagram. In this manner one generates 
a single two-loop diagram with a symmetry factor in front of it: One can then 
compare this to what one would get from a standard Feynman diagram calculation. 
Similar formulae hold, of course, for the crossed versions of these diagrams. 

After checking that the method is correct, the second question to be answered is 
how one can do Feynman loop integrals with a nonlocal vertex in it. In particular, 
how does one perform the following integral: 



{27lY (M2 - /2)(M2 - (pi + p2 - ■ 



(2.49) 



First we remark that this integral has to be calculated anyway, even if one would 
try to calculate directly the two-loop diagrams. It corresponds to the fish diagram 
in Fig. where the divergence of the one-loop subdiagram has already been 



4|c5. As can be seen here, one 
le calculation at a stage where 



cancelled by the insertion of the counterterm Fig. 
of the advantages of the method is that one starts tl 

the one-loop subdivergence has already been subtracted.^ This fact is of relevance 
in particular in CHPT, where one may thus first do the renormalization of the off- 
shell one-loop amplitude before diving into the forest of two-loop integrals. Next we 

^The subtraction of the subdivergence in the diagram Fig. ^3 is much easier to accomphsh, 
because the corresponding integral is a product of two one-loop mtegrals. 
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note that the loop function J(s) admits a dispersive representation in d dimensions 
m, 11: 



J{t) = t r -fe^ , (2.50) 

^ ^ JiAP X{x-t) ^ ^ 

where the measure [dx] is given in appendix ^ (for = 1). By inserting this 
expression into Eq. ( ^.491 ) we see that the problem is transformed into the calculation 
of a one-loop integral with three propagators, where one of the propagators has a 
variable mass over which we finally have to integrate: 

[dx] f dH {l-p,f ^2.51) 



4M2 X J (27r)^(M2-/2)(M2-(/-p3-p4)2)(x-(/-J93)2) 

We discuss in appendix ^ how this integral can be done explicitly . 



Finally, we are left with the graphs that we have not yet taken into account: The 
tadpole and counterterm insertions displayed in Fig. |5di,...,4 can easily be taken care 
of by simply replacing in the one-loop integrals the parameter M with the physical 
mass Mp. The contributions from wave function renormalization is taken care of by 
multiplying the tree graph with 

In summary, the advantages of the approach just described are that i) by using 
the renormalized one-loop amplitude as a single nonlocal vertex, the divergences due 
to one-loop subgraphs in nonfactorizable integrals never appear; ii) the diagrams 
that produce the mass renormalization inside the one-loop amplitude are taken into 
account trivially and do not have to be calculated explicitly; iii) the method also 
applies to the calculation of form factors, or any process where one part is given 
by an elastic two-body process like 77 — > vrvr [Q; iv) the method can easily be 



implemented into a computer algebra program like FORM |25| 



3 CHPT to two loops 

In the following we discuss the methods used in Ref. [|10| to evaluate the elastic tttt 
scattering amplitude to two loops in the framework of CHPT. This section contains 
details concerning renormalization, whereas the two-loop expressions for the pion 
mass, the pion decay constant and the scattering amplitude are given in the following 
section. 

3.1 The lagrangian 

The effective lagrangian consists of a string of terms, 

= £2 + + h'^C.Q + .... (3.1) 

Depending on the application one has in mind, one has to keep external fields in C^g. 
In the present case, we wish to calculate Mt^^F.,^ and the vrvr scattering amplitude. 
The pion mass and the scattering amplitude may both be calculated from on-shell 
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quantities, and external fields are thus not necessary. On the other hand, the pion 
decay constant needs a weak current as an external probe - it is an off-shell quantity, 
not accessible through on-shell matrix elements. The use of the axial current to 
evaluate to two loops is cumbersome because of the presence of Lorentz indices. 
In the following, we use instead the two-point function of pseudoscalar densities. A 
Ward identity relates the residue of this quantity to the pion decay constant. 

There are many choices for the pion fields to be used in the effective lagrangians. 
Of course, the result for physical quantities is always the same. A convenient choice 
to minimize the number of diagrams is the sigma model parametrization. In the 
following we work in the isospin symmetry limit m„ = m^. We have 



C2 = ^{uX + X+) , (3.2) 



with 



U = a 





^2 , 








V2n- 





X+ = u^xu^ + ux^u , 
X = 2B{ml + ip) , m = -{niu + nid) , (3.3) 

with = U. Furthermore, p = p^r'^ is the external pseudoscalar field referred 
to above. The symbol {A) denotes the trace of the two-by-two matrix A. The 
lagrangian £4 is 

£4 = jZkP^ + ■■■ , (3.4) 



i=l 



where 



(3.5) 



with 



X- 



u^d^XU^ - ud^x^u . (3.6) 



The ellipsis in Eq. ( p.4|) denotes terms that do not contribute to the quantities con- 
sidered here. The low-energy constants li are divergent and remove the ultraviolet 
divergences generated by one-loop graphs from £2. We discuss them in more detail 
below. 

The complete effective lagrangian is not yet available: Whereas the list of the 
necessary counterterms has been published p6| , their divergence structure at d = 4 is 
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still under investigation p^. The knowledge of these would provide us with a check 



on the calculation, because these divergences must cancel the ones that we find in 
the two-loop calculation. Apart from this check, however, that analysis would not 
be of further use in the present context, because the scale-dependent finite pieces of 
those couplings are still largely unknown. Nevertheless, the structure displayed in 
26[] shows that there are no algebraic constraints between the counterterms at 0{p^^ 



needed in the expressions for the pion mass, pion decay constant or tttt scattering 
amplitude. 

For the following considerations, it is useful to recall that the expansion in powers 
of the momenta is equivalent to an expansion in inverse powers of F^. In an obvious 
notation, the chiral expansion for the pion mass, the pion decay constant and the 
elastic tttt scattering amplitude is 

Ms An) = + + + (3.7) 

The one-loop contributions m4, /4 and A4 have been determined in - here, we 
calculate me, /e and Aq. Therefore, we need to keep in the effective lagrangian only 
contributions up to 0{F~^). We write symbolically 

A = >Cfci„ + ai— + aa-^ + 03— + 0(F ) , 

A = ^^Ji + ^^^ + 0{F-') , (3.8) 

where the dimension of the couplings is 

[oj] = mass^ , [bi] = mass^ , [cj] = mass^ . (3.9) 
In order to generate a term of order _p-2m^ ^^^^ consider all products 

^ni'Ca'Cs' ; ^1^1 + "^2^2 + "^3^^3 = m . (3.10) 

Furthermore, it is very convenient to collect the terms quadratic in the pion fields in 
the kinetic part and to expand in powers of bi and C2 only afterwards. As a result, 
one has to consider the following products: 

me, /e : aj; ag; 62 

Ad : af; 0102; as; 0162; b^; C3 . 
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3.2 Renormalization at two— loop order 



3.2.1 General analysis 

The loop contributions to Mt,,Ft, and to the scattering amplitude A{s,t,u) are 
divergent in the limit d —>■ A. Using the notation ( |3.7|) , we are concerned here 
with the renormalization of the quantities tuq, Jq and Aq. As is guaranteed by 
general theorems of renormalization theory |2^, the divergent parts of tuq and /e 
are mass independent in dimensional regularization. Likewise, the divergent part of 
Aq is a polynomial in the external momenta and in M^. The most general crossing 
symmetric polynomial arising at order is 

^6^po, ^ A^^)M'^ + A^'\sM^ + A^''h^M^ + A^^\t-ufM^ 
+A^^h' + A^'^hit - uf 

with six coefficients A^^\ . . . , A^^\ Like itiq and /e, they receive contributions from 
two-loop graphs with C2, one-loop graphs with one vertex from £4 and tree graphs 
generated by C2 + ^^4^ + Cq, see Eq. ( p.ll| ). The following analysis applies to each 



of the eight coefficients mg, /e, A^^\ . . . , A^^^ separately. Denoting these coefficients 
generically as Q, we write 

Q — Qioop ~l" Qticc ■ (3.12) 

Since Q is part of a measurable quantity, it must of course be finite (and scale 
independent). The two components Qioop and Qtrco, on the other hand, are both 
divergent. Concentrating first on Qioop, dimensional regularization produces for 
and /e the general form 

Qioop = Tm i^fMx{d) -J2li{d)yi{d)^ , (3.13) 

where x{d),yi{d) are dimensionless functions of d, finite at d = 4. (For the sake of 
clarity, we exhibit in this section the dependence of the low-energy couplings li{d) 
on the space-time dimension d.) The same structure is found for the six coefficients 
^(*) We perform a Taylor series expansion, 

x{d) = Xq + XiW + X2W^ + 0{w^) , 

Viid) = yio + yiiw + yi2w'^ + 0{w^) ;i = l,...,4; w = ^-2 



with real numbers Xo,Xi,X2,yio,yii,yi2- The tadpole integral Tm (|2.17|) and the 



low-energy constants li{d) are both diverg ent as w — * 0. The Laurent expansion of 
the tadpole integral is displayed in (|2.18|) . For the low-energy constants li{d), we 
write 

k{d) 



(4vr; 



2 



12 1 
7i = 3 ' ^2 = - , 73 = "2 > 74 = 2 . (3.14) 



19 



For reasons already explained in Sect. 2 in the framework of the A^-component 0^ 
theory, we have not expanded the renormalized couplings around w = 0. 

See comments after Eq. ( p.24| ) and the discussion below. 

Since there must not be any terms of the form 

InM/// 

w 

in Qioop, there are eight consistency conditions (for each of the eight coefficients 
m6,/6,A«) of the type i 



Xo 



4^ 



iO 



(3.15) 



i=l 



This absence of mass singularities in the divergences provides an extra check on our 
calculation. Using the equality ( p.l5|) that relates one- and two-loop parameters, 
one gets (from now on, the summation convention is implied for i = 1, . . . , 4) 



Qloop 

Q2 

Qo 



i ^ H h Qo + 0{w) \ 



(47r)4 

1 

-Xo , Qi = xi- lf^{^i, w)yio - -liVii 
xoa{M/i^y + 

+X2 - /,>(/i)yil - -7il/i2 



2x,-l^^{fi)y.o 
1 



a{M/fi) 



where 



/-(/i) = /-(/x,0). 



It is seen that the coefficient c(M/ fi) in the Laurent expansion ( p.l8|) of the tadpole 
integral Tm does not enter due to the relation ( p.l5|) . 

The pole terms can now be absorbed by the tree-level contribution Qtroc which is 
a certain combination (depending on the specific observable under consideration) of 
coupling constants in the general chiral lagrang ian C q. Denoting this combination 
generically as z{d), the expansion analogous to (|3.14| ) is 



QtrUd) = z{d) 



Aw 



{-% - — + (4vr)\-^(/.) + 0{w) \ . (3.16) 



{Any 



The sum Qioop + Qtrcc is finite and independent of the scale by construction. Taking 
the limit 4 yields 



lim Q 



(4vr) 



(3.17) 
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Recalling the scale independence of Q,TM,li{d) and of z{d), one derives the renor- 
malization group equation for the renormalized coupling constant z^^: 



dfi 



(4vr) 



2xi - lf^{^i)y^o - -iiV, 



il 



(3.18) 



Comparing Eqs. ( p.l9|) , ( p.l5|) and (|3.17|) , we recover the well-known fact [||, |2^, 
|T9| that the coefficient of the double logs can be determined solely from one-loop 
diagrams with a single vertex from £4. 



3.2.2 Modified minimal subtraction 



The above renormalization procedure corresponds to m inim al subtraction, where 
powers of In 47r + r'(l) occur in the final expressions, see (|2.19|) . These terms can be 
absorbed by splitting from the renormalized couplings If^ilA ^-nd z^^^n) appropriate 
finite pieces P3. The procedure is based on the relation 



exp w[V\l) + ln4vr] = exp £ 

V^'^) n=2 ^ 



(3.19) 



which shows that a simple factor is responsible for these terms. To remove them, 
one pulls out a factor c^^ in the definition of the renormalized couplings in (|3.14| ), 



h{d) 



(/ic) 



2w 



(47r)2 {2w 



(3.20) 



as a result of which Qioop in (|3.13|) becomes 



Qioop = if^c)^"" { {yic)-^'"flx{d) - (/ic) 



-2w 



4 r 



E 



Pulling out in the analogous manner c^'^ in Eq. ( 3.16 ) 

(/ic)^"' 



l/,(c/) .(3.21) 



(3.22) 



shows that these redefinitions amount to the change 



In 



M 



In In c 

/i 



(3.23) 



in Eq. (|3T7D . 

The traditional choice for c in CHPT is 

1 



lnc= — [ln47r + r'(l) + l] . 
2 
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In this scheme, one uses the notation 



^^(d) = + - ^[l^^vr + r'(l) + 1]} + 0( 



where 



(47r)2 

according to ( |3.20|) . Similarly, we write 

z^ifi) = (3.24) 

for the above choice of the constant c. 

We have thus arrived at the final expressions that will be used for M^, and 
A{s, t, u): 



Q 



(4vr) 







2 








+ 













\2ir/ 



'MO 



l + 21n 



+X2 - -liVii - (47r)2/[(/x)?/a + (47r)S"(/i) 




The scale dependence of the renormalized couplings is 



dfj, 



(47r)4 

7£_ 

(4vr) 



[2xi - (An) ri{fi)yio - liVn] , 



3.2.3 EOM and Laurent expansion of the coupling constants k 

We include here a brief discussion of some technical aspects of the renormalization 
program at 0{p^). In Sect. 2, it was shown that the counterterms of 0^ theory can be 
written in different forms using the EOM. Transforming from one set of counterterms 
to another leaves the functional Zi unchanged and produces an additional local 
action of 0{h^) that can always be absorbed by changing the coefficients of the 
counterterms in Z2. 

The situation in CHPT is similar. Using the EOM of CHPT, both the lagrangians 
of 0{p^) and O(p^) can take different forms. Since toO(/) the la grangian Cq enters 
only at the classical solution, such a modification has no effect at all as far as Cq is 
concerned. The same is true for all quantities of O(p^). However, the couplings in 
the lagrangian £4 also appear in one-loop and tree diagrams contributing at Oijf'). 
In this case, different forms of £4 lead to different results a priori. 

As for the 0^ theory in Sect. 2, the natural place to discuss those modifications 
is the generating functional of Green functions. Referring to a forthcoming paper 
27[| for an explicit proof, we present here only the final result. Lagrangians of 0{j)^) 
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that differ by the EOM lead to generating functionals that differ by local actions 
of 0{p^). Since those actions are of course chirally symmetric they can always be 
absorbed by a redefinition of the low-energy constants of 0{p^). 

If the low-energy constants are determined phenomenologically by comparison 
with experiment, the result is of course independent of the form of £4. On the other 
hand, using a model for the constants of 0{p^) like resonance saturation does not 
specify the form of £4. Therefore, different forms of £4 lead to different numerical 
results because the constants of 0{p^) are by definition unchanged. For the vrvr 
scattering amplitude, it turns out that the coefficients 61,. . . ,64 are affected by this 
ambiguity while 65 and are not modified. We come back in Sect. 5 to a numerical 
discussion of this ambiguity. 

The second issue we want to address is the dependence of the low-energy con- 
stants li of 0{p^) on the dimension of space-time. In Eq. ( |3.14| ), we have defined 
functions lf^{fi,w) that contain all non-singular pieces of the li{d). These functions 
enter in the quantit y Qi in Qioop and they are subtracted by the appropriate decom- 
position of Qtrcc in ( |3.16|) . As a consequence, only the coefficients = /f^^(/i,0) 



appear in renormalized quantities. 

Although this is a completely legal subtraction procedure, it is by no means 
obvious that this procedure can be applied consistently for all possible processes. 
In other words, the question is whether the low-energy constants of 0{p^) can be 
decomposed as in Eq. ( p.l6| ) in a process-independent way such that only the IfrM 
appear in all observable quantities. To discuss this issue for the case of minimal 
subtraction, we expand the li{d) up to first order in w: 



2w 



(4vr) 



^ + lT^{fi) + S,{fi)w + 0{w^ 
2w 



(3.26) 



The terms of O(iy^) cannot contribute to 0{p^) because the k appear only in dia- 
grams that are not more singular than 1/w (except for the divergent parts of the 
li themselves). However, the coefficients 5i(/i) certainly do appear at 0{p^). The 
subtraction procedure ( p. 161 ) leads to the following two questions: 

• Do the coefficients 5j(/i) appear always in local contributions of 0{p^) only? 

• If yes, can they be absorbed into the coupling constants of 0{p^) in a process- 
independent manner? 



The crucial question is of course the first one. It can be shown [^T], ^ that 
the terms involving the 6i indeed take the form of local actions in the generating 
functional of 0{p^), see also the analogous discussion in the framework of the A^- 
component 0^ theory in Sect. 2. Although these contributions depend on the choice 
of the lagrangian of 0{p^) (EOM), they can therefore always be absorbed in the 
low-energy constants of 0{p ) in a process-independent way. Thus, the subtraction 
procedure (|3.16| ) is a special case of a general method. It has already been used in 



4 Elastic TTTT scattering to two loops 
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Figure 6: Butterfly diagram contribution to the two-point function ( ^.11 ). 



4.1 Pion mass 

We consider tlie connected two-point function 



d xe 



-ip{x-y) 



(O|T0Xx)0'=(y)|O), 



(4.1) 



at two-loop order. According to the formula (|3.11| ), the topology of the graphs is the 
same as in 0{N) theory, except for the term proportional to (j)^. We may therefore 
again simplify the calculation by evaluating the tree and tadpole diagrams displayed 
in Fig. ^ with an appropriate choice for the propagators and for the interaction 
vertex, and to add the sunset diagram Fig. ^lo, using the vertices that correspond 
to £2- Finally, the term from 02 generates the butterfly diagram Fig. ^ Collecting 
these contributions, we find 



Ml - p2 



ik 



(4.2) 



where 



Mi 



M' + 



M^ 



F2 



2/3 + 



M 



2M2 



1 - ^ + ^pT^A/ 



M^ 



+ 0{M^ 



Q^Tm 



i=l 



rz 



+ 0{M'' 



(4.3) 



where p4| 



and 24 



Qf 



- (204 - 632u; + 1263^;^ + 0( 
96 J 

^ [28 - 30w + 31w^} + 0{w^) , 

8 - lOw + llw^ + 0{w^) , 
3 - 4w + 4w^ + 0{w^) , 



(4.4) 



— (96 - 464W + 1185wH + 0{ 
96 J 

Qi ) Q2 ~ Q2 ) Q'c 



w 



2 . 



(4.5) 
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The function R{p^) in the propagator is regular at = M^. It is not needed in the 
following, and we no not display its explicit form here. The quantities tm and 
rz denote the counterterm contributions from the lagrangian Cq. We note that tm 
renders the pion mass finite at two-loop order. 



4.2 Pion decay constant 

We consider the correlators of two pseudoscalar currents 0, 
q = {u,d) 



(4.6) 



at two-loop order. Apart from the overall normalization, G^^ differs from the two- 
point function A*'^ only by terms generated by the low-energy constants I3 and U 
(and, of course, by the counterterms from Cq). We find 



Gl 



Ml - p^ 



ik 



(4.7) 



where 



Gl 



G^ 



M2 , 
1 + — ^ 2/4 - 4/3 



F2 



^Grh 



i=l 



+ 



{ll + 8/4/3 + 4/3' + re} 



+ 0(M^ 



(4i 



with 



G = 2FB , 

Q — Q 1 Qi — Qi ) Q2 ~ Q2 5 

= 2{8-7w + 7w^) + 0{w^) , 

= l-w + w^ + 0{w'') . 



(4.9) 



The counterterm contribution from Cq is denoted by r^. It renders finite at two- 
loop order. The regular part Rj^ is not needed in the following, and we therefore do 
not display it here. 

The relation between the divergence of the axial current and the pseudoscalar 
density implies 



F^M^ = mG, 



(4.10) 



The above expressions for the pion mass and for the residue G^^ therefore allow us 
to calculate the pion decay constant to two loops: 



1 + 



F2 



Tm 
F2 



71 



i=l 



(4.11) 
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Figure 7: Contributions to elastic vrvr scattering in CHPT. The graphs displayed have 
a different topology than the ones occurring in 0^ theory. The filled box denotes a 
contribution proportional to the low-energy constant I3. Graph ai is included by 
use of the formula indicated in Fig. ^, whereas graphs a2 ■ ■ ■ as are not. These must 
be taken into account separately. 



with 

= — — [240 - 656w + 1125w'^) + 0{w^) , 
192 

Qf = -^Q? , Q2 = -\Q2 , Qa^ = 2 , = hQ^i . (4.12) 

The counterterm contribution from is denoted hj rp. It renders finite at 
two-loop order. 



4.3 The TTTT scattering amplitude 

We have discussed in Sect. 2 how one can simplify the calculation of the elastic 
scattering amplitude in the iV-component 0^ theory by using the renormalized one- 
loop amplitude as a single nonlocal vertex. In case of CHPT, there are a few 
modifications to be taken into account: 



As we discussed already in the case of the pion mass and pion decay constant, 
there are additional diagrams to be taken into account, see the formula ( p.ll| ). 
Some are displayed in Fig. ^ The first of these is taken care of by using the 
formula illustrated in Fig. ]g, whereas the tadpole (Fig. [^2), the butterfly 
(Figs. ^3,a4) and sunset diagram (Fig. [^as) have to be added by hand. 

Since the effective lagrangian in CHPT has derivative interactions, the vertices 
entering the calculation have a nontrivial off-shell behaviour. This means that 
the one-loop amplitude entering the above formulae must first be calculated 
with the appropriate off-shell legs in d dimensions. Since the relevant expres- 
sion is not available in the literature, we display it in appendix 0. 



Needless to say that, despite the labour saving organization of the calculation that 
we just described, it is still rather long and tedious. 

In analogy to what we have presented above, we might now display the amplitude 
in d dimensions. Since the expression is rather long, we prefer to display only the 
renormalized final expression. 
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We use the notation 

(vr"(p4)7r^(p3) ont\7i%p^)n\p2) in) = {Ap,)7r%p3) m\n^{p,)7r'{p2) in) 
+i{27rys\Pf - Pi) {6''^6^'^A{s,t,u) + cycl.} , 

where s, t, u are the usual Mandelstam variables, expressed in units of the physical 
pion mass squared M|, 

s = (Pi + P2f/Ml ,t^{ps- p,f/Ml ,u^{p,- pif/Ml . (4.13) 

Using these dimensionless quantities, the momentum expansion of the amplitude 
amounts to a Taylor series in 



X2 = 



Ml 
F3 



where F^^ denotes the physical pion decay constant. We find 
A{s, t, u) — X2[s — 1] 



-X2 



+ X2' 
+X2' 



h + b2S + hs^ + h^{t - uf'\ 

FW(s)+G'«(s,t) + GW(5,w) 
hs^ + hs{t - uf'\ 



(4.14) 



with 

F«(s) 

F^'\s) 



1) , 



^J(t)(14 - As 



lot 



St + 2r 



1 /503 o 



929 



S87 _ 140 

27^~ ~9~y 



G^'\s,t) 



"^^ * I 167rm08 ' 54 
+ 61 (4s -3) + b2{s^ + As - A) 

+ ^ ( 8^2-21 + 48s- 32) + ^ (168^-71^2 + 112s-48)j 

+ ^Ki{s) |20s^- 119s^ + 210s- 135- :^7r^(s- 4)1 
18 I 16 J 

+ li^2(s) {s7r2-24} + ^X3(s) {3s2-17s + 9} , 

/267 



32 

J{t) 



167r2 



412 s 
^~54 



(r + 5t+ 159)-t 



V216 



„ 727 1571 \ 
t - —t + 



108 



108 



+ 6i(2-i) + ^(i-4)(2i + s-5) - -(i-4)2(3i + 2s-8) 
3 6 



27 



+ ^ (2s(3t - 4) (t - 4) - 32t + AOt^ - llt^ ) | 

1 f TT^ 

+ —Ki{t)\ 174 + 8s-10t3 + 72t2-i85t (t-4)(3s-^ 

36 [ 16 

+ ^K^it) |l+4s+^t(35-8) I 

- Ksit) {l + 3st - s + St"^ - 9t} + - lUit) {A - 2 s - t} . (4.15) 
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The loop functions J and Ki are 
/ J \ / 



z 

z 


\ Nzs-^ 7r2(iVs)-i 



-AN \ 



TX 



( \ 

h 



and 



^ 2''^ + 3^^+iV^ + H2iV^ 



where 



In 



i-1 



A^Vi 1 



1 - 



= 167r2 



The functions s~^J and s^^Kj are analytic in the complex s-plane (cut along the 
positive real axis for s > 4), and they vanish as |s| tends to infinity. Their real and 
imaginary parts are continuous at s = 4. The coefficients hi in the polynomial part 
are given in appendix 0. 



Finally, we compare our result with the calculation performed in Ref. 0. There 
are two basic differences between the two works. First, the present calculation is 
done in conventional chiral perturbation theory p[ , whereas the authors of work 
in a scenario in which the quark condensate may be small or zero |^ (generalized 
chiral perturbation theory). The second difference concerns the use of a lagrangian 
framework in the present approach, whereas S-matrix methods are applied in 0: 
Starting from the expression for the amplitude at O(p^), unitarity allows one to 
determine the absorptive part at O(p^). One can then construct an amplitude that 
has exactly these absorptive parts. Of course, this procedure does not determine 
the polynomial part at this order. To illustrate, the diagrams of Fig. |^2,3,4,5 do not 
have an absorptive part and can therefore not be determined in this manner. There 
are additional graphs of this kind. In the amplitude given in 0, these diagrams - 
and further polynomial contributions from graphs that do have an absorptive part 
- are lumped into a polynomial 



uY + a^s^ + aQs{t 



(4.16) 
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with a priori unknown coefficients that depend on the pion mass and on the low- 
energy constants. 

With this procedure, one does not have to work out individual Feynman dia- 
grams, because only the total absorptive part at OijP) - fixed through unitarity by 
the amplitude at O(p^) - is needed. In addition, the expansion of the pion mass, 
of the pion decay constant and of the wave function renormalization need not be 
worked out at 0{p^). 

On the other hand, our result ( |4.14| , |4.15| ) at 0{p^) - obtained by painstakingly 
evaluating in the manner described above all the Feynman diagrams generated by the 
effective lagrangian ( |3.1|) - reveals in addition to the findings of the dependence 
of the six coefficients ( |4.16| ) on the pion mass and on the low-energy constants of 
both 0{p'^) and O(p^), see appendix 0. In particular, we find that these coefficients 
contain mass singularities that are known to be important numerically both at O(p^) 
and at 0{p^) [|l^, |T0| (see also Refs. |2^, |4|, |33|). Furthermore, the knowledge 



of the mass dependence of the amplitude allows one to evaluate all quantities even 
at unphysical values of the quark mass, and to confront the result with lattice 
calculations |0. This is not possible with the amplitude provided in Ref. 0. 

We have checked that the amplitude of (restricted to standard CHPT) and 
the field theoretic calculation presented here agree as far as the absorptive part of 
the amplitude and the general structure of the real part is concerned. 



5 Numerical analysis 

The elastic tttt scattering amplitude is expressed in terms of the external momenta, 
the physical pion mass, the physical pion decay constant, and the coefficients 6j, 

A{s,t,u) = f{pi, . . . ,p4; M^, F^;bi, . . . ,be) . 

Before one can perform a numerical analysis, one thus needs an estimate of the 
low-energy couplings bi. The key relations ( p.l|) used in [|10| for that purpose fix 
the bi in terms of 

1 

- chiral logarithms L = In — ^ , 

- the low-energy couplings IKfi), . . . , lUfi) from £4 , 

- the low-energy couplings r[(yu), . . . ,rg(/i) generated by Cq . 

For a given set of the low-energy constants Z[ it, therefore, suffices to estimate 
the new couplings r[ at 0{p^). To achieve an order-of-magnitude estimate, we 
use a method that has been successfully tested at O(p^) 0, |^: We work out the 
contributions from exchanges of heavy states to the scattering amplitude and assume 
that these effects account for the bulk part in the low-energy couplings at 0{p^). In 
particular, we include the effect of meson resonance exchange with masses smaller 
than 1 GeV, i.e., vector and scalar resonances. In addition, we also consider kaon 
and eta contributions of 0{p^). We describe in the following subsection the relevant 
resonance couplings and work out the corresponding values for the r^. 
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5.1 Resonance saturation 



5.1.1 Vector meson exchange 

At O(p^), the contributions of vector meson resonances to the low-energy constants 
are obtained most naturally with the tensor field representation of spin-1 mesons 
0, The situation is different at 0{p^) where the more common vector field for- 
malism produces directly the relevant counterterm couplings, as observed previously 
for Compton scattering on pions and for 77 — tttt Q . We follow the same pro- 
cedure here and comment later on the difference to the tensor field representation. 

The relevant couplings of the vector field V^, representing the octet of vector 
mesons (the singlet does not contribute to vrvr scattering) , to the pseudoscalar mesons 
are given by [061 



= (v;.K,^i) + {v,K,x-]) (5.1) 



with real dimensionless coupling constants gv,fx- switch to chiral SU{3) here 
(and also in the next subsection for the discussion of scalar exchange) to determine 
both couplings gy, fx from vector meson decays. From the Lagrangian ( |5.1| ) one 
finds that the combination 



V 



determines the amplitude for the decay of a vector meson into two pseudoscalars. 
From the experimental widths for p — tttt and K* Kn one obtains Qvfx < 
with 

gv = 0.09 , = -0.03 , (5.2) 

if we define gv to be positive. V exchange on the basis of ( [5.1[ ) gives rise to three 
types of contributions proportional to gy, gvfx /x' respectively. With 

using My = Mp for the numerical values, we obtain the following V exchange 
contributions to the low-energy constants rj: 

r\ = — 16ay — I6by — Acy 
= 20ay + 16by + 3cv 

= —lay — 36y (5.3) 



-2.2 ■ 10" 
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ay + hv 
3 



In the tensor field representation, there is a single coupling (relevant for tttt scat- 
tering) between vector mesons and pseudoscalars at lowest order p|, |3^. Up to 
normalization, the corresponding coupling constant is given by gy ■ For instance. 



expanding the p exchange amplitude of Ref. p[ to 0{p ) is equivalent to (|5.3| ) with 
6y = cy = 0. Higher-order couplings can only contribute at O(p^) through inter- 
ference with the gv amplitude. Thus, the contributions involving that appear 
naturally in the vector field formalism would have to be added by hand in the tensor 
field representation as explicit local counterterms. 



5.1.2 Scalar meson exchange 

To estimate the effect of scalar resonances for the r^, we take the lowest-order 
lagrangian of Ref. for octet and singlet scalar fields 5, 5*1, 

^5 = Cd {Su^'u^) + c„ {Sx+) + ^5*1 (n^'M^) + c;;^! , (5.4) 

and expand the resonance exchange amplitudes up to Oijf'). The resulting scalar 
contributions to the r,- are: 



rf = 





rf = 


3^aCm Ca) 1 [ 


rf = 


3Mr"^^"' ' Ml 


rf = 





rf = 




rf = 


. 



Cd 



(5.5) 



Note that the physical pion mass and the physical pion decay constant receive con- 
tributions from scalar exchange as well. 

For the numerical evaluation we use the same values for masses and coupling 
constants as in Ref. [03: 



Ms = Ms, = 983 MeV 
= 42 MeV , Q = 32 MeV 

Cin = Cm/ VS , Cd = Cd/VS . 
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5.1.3 Kaon and eta contributions 



In the framework of chiral SU{3), K and 1] mesons contribute at O(p^) to tttt scat- 
tering via loop diagrams. Restricting the scattering amphtude of 0{p'^) evaluated in 
chiral 5*^7(3) to SU{2) x SU (2) by an expansion in inverse powers of the strange 
quark mass, one arrives at the following contributions to the rj due to K and t] 
mesons: 



K 



K 



^5 



31F! 



57607r2Mi 



23047r2M|. 



29F! 



76807r2M|. 



25607r2M|; 



(5.6) 



23FJ 



153607r2M|: 

TT 

153607r2M|^ 



5.1.4 Numerical values of the 

Putting the various contributions of the previous subsections together, we obtain 
the following numerical estimates for what we shall call the resonance contributions 
to the low-energy constants of 0{p^): 

rf = -0.6 ■ 10"^ 

= 1.3 -10-" 
rf = -1.7-10-^ 

rf = -1.0 -10-^ (5.7) 

= 1.1-10-^ 
rf = 0.3 ■ 10"^ . 

These values are the ones used in Ref. [|lOl. They indicate an 0{p^) version of vector 
meson dominance: Despite some strong cancellations (especially for , r^), all six 
constants rf are dominated by the exchange of vector mesons.^ 

Of course, the constants 6j become scale dependent if we now replace the r[(/i) 
by the above rf . At ©(p"^), the empirically observed resonance saturation of low- 
energy constants amounts to the statement that the resonance contributions 
describe the phenomenologically determined coupling constants quite well for a scale 
fi between 500 MeV and 1 GeV. Pending a more refined analysis along the lines 

^ Other estimates of the low-energy constants of 0{p^) can be found in Refs. pq]. 
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discussed below, we assume for the time being that the same approximation is 
meaningful also at 0{p^). Of course, only observables that are relatively insensitive 
to scale changes in the range mentioned can be estimated in a reasonable manner 
in this way. 



5.2 The constants 6, 



Together with the couplings /[ at 0{p^), the estimates (^.7|) allow one to work out 
the constants &i, . . . , b^. We used in Ref. |T0| values for the /[ that were determined 
mainly from an O(p^) analysis 0, with additional input from a dispersive estimate 
of higher-order effects in Ku decays . As emphasized in our previous paper, all 
those determinations were faced with the problem that the /j are mass independent, 
whereas the physical quantities, from which the /J were determined, include quark 
mass effects. We will illustrate below the relevance of those quark mass effects by 
example. 

With the scale-independent couplings k (cf. App. D) found in p 



h 
h 



-1.7 
2.9 



6.1 
4.3 



(5. 



and with the constants rf^ in ( |5.7| ), one arrives^ at the bi displayed in table |l|. 
There, we have split the contributions to the bi into one-loop and two-loop effects. 
Furthermore, we indicate separately the contributions from the rp. The values 
shown correspond to the scale /i = 1 GeV [fi = 500 MeV in brackets). (We display 
the values to two digits for ease of reproduction of the numbers.) Below we refer to 
the bi in this table as set I. The following remarks are in order. 



Table 1: The constants bi (set I). We use Eqs. ( ^.7| , |5.8|) as input to evaluate bi 
from Eq. ( p.lj ). The numbers correspond to /i = 1 GeV, in brackets we display the 
corresponding values at = 500 MeV. 





1-loop 


2-loops, ri=0 


from rf- 


total 




-7.34 


-1.76 (-1.82) 


-0.01 


-9.11 (-9.17) 


10% 


6.74 


2.07 ( 2.12) 


0.03 


8.84 ( 8.89) 


10%3 


-0.84 


-3.05 (-3.43) 


-0.38 


-4.27 (-4.65) 


10%^ 


5.56 


1.72 ( 1.71) 


-0.22 


7.05 ( 7.04) 


10% 




1.22 ( 1.24) 


1.10 


2.32 ( 2.34) 


10% 




1.19 ( 0.77) 


0.30 


1.49 ( 1.07) 



As in Ref. |10 , we use = 93.2 MeV, = 139.57 MeV. 
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The change induced by /i = 1 GeV ^ jj, = 500 MeV is of the same order 
or bigger than the contributions from rf (see also table with the possible 
exception of 64. 

Besides the scale dependence, there is yet another source of uncertainty that 
has to do with using the EOM in the lagrangian of O(p^). As discussed in 
Sects. 2 (for 0^) and 3 (for CHPT), this ambiguity can always be resolved 
by a redefinition of the coupling constants of O(p^) - physical observables 
are of course unaffected. Although this is a completely general result for the 
generating functional to O(p^) we nevertheless face a practical problem 
here because our estimates of the r[ via resonance exchange are not only scale 
independent but they also make no reference to the form of the lagrangian of 
O(p^). We have checked that using the O(p^) lagrangians of either Ref. 0] or 



of Ref. 1^ induces a change in 2,3,4 by an amount that is smaller than the 
change due to a different choice of scale. The constants 65^6 are unaffected by 
this procedure. 

The main uncertainties in the hi stem from the couplings of O(p^), as we now 
illustrate. In table |, we have displayed the contributions from U to hi at one- 
loop order. They are at least an order of magnitude bigger than the ones from 
the - therefore, uncertainties in the latter are swamped by uncertainties 
in the /j. One of these are quark mass effects^ rnentioned above. To illustrate 
these, we repeat the analysis of Ref. [0] to fix /i, I2 from the D-wave scattering 
lengths a^, but now to 0{p^). Keeping [3, fixed, the experimental values 
given in table § lead to 

h =-1.5 , , 

h = 4.5 

if I3, 1 4^ from ( |5.8| ) are used as input. The hi that result from this exercise are 
displayed in table || and referred to as set II in the following. Since 61,2 do not 
depend on I21 they remain largely u naff ected, whereas the change in 63^... ,6 is 
seen to be substantial. The values ( ^.9|) should be compared with the values 
obtained at 0{p*), li = —2.3 and I2 = 6.0, from the same input [@]. Values 



similar to (p.9| ) were found in |2^ comparing a Roy equation fit to our p 
amplitude ||lO|| . 

Similar effects are expected from changes in I4 that does contribute substan- 
tially to 61^2 • As we already mentioned in our previous work a more 



detailed analysis is therefore needed to obtain reliable values for the li 
tent with an analysis at 0{p^), and, at the same time, values for the hi with 
error bars attached. Such an analysis is under way |^ . 



This brings us to the work of Ref. |T3] who find from dispersion sum rules the 
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Table 2: Contributions from the individual U in Eq. ( [5 .81 ) to the constants bi in the 
one-loop approximation. Note that 65,6 start at two-loop order. 



1-loop 


k 


k 


^3 


k 


analytic 


total 


10% 


-1.44 





-0.92 


-5.44 


0.46 


-7.34 


10% 


1.44 








5.44 


-0.14 


6.74 


10% 


-3.59 


6.44 








-3.69 


-0.84 


10% 





6.44 








-0.88 


5.56 



Table 3: The constants 6j (set II). We use Eqs. (|5.7|J5.9| ) as input to evaluate bi 
from Eq. ( U.l ). The numbers correspond to /i = 1 GeV, in brackets we display the 
corresponding values at = 500 MeV. 





1-loop 


2-loops, rj=0 


from rf^ 


total 


10% 


-7.17 


-1.44 (-1.62) 


-0.01 


-8.63 (-8.80) 


10% 


6.57 


1.41 ( 1.71) 


0.03 


8.01 ( 8.31) 


10%3 


-2.11 


-0.14 (-1.85) 


-0.38 


-2.63 (-4.34) 


10% 


3.87 


1.17 ( 1.20) 


-0.22 


4.82 ( 4.85) 


10% 




-1.47 (-0.22) 


1.10 


-0.37 ( 0.88) 


10% 




0.53 ( 0.42) 


0.30 


0.83 ( 0.72) 



values 



be 



-3.7 ±2.4) ■ 10-3 
(4.8 ±0.3) • 10-3 
(1.4 ±0.6) ■ 10-^ 
(1.0 ±0.1) ■ 10-^ 



(5.10) 



Within the error bars quoted, these values are consistent with our set II. The 
theoretical framework used in [0 does not allow these authors to determine 
bi^2 ^ rather, these quantities can in principle be determined from precise data 
on elastic tttt scattering, allowing then a determination of the size of the quark 
condensate in the chiral limit. 

vi) Recently, Wanders has also extracted the 63,4,5,6 using crossing-symmetric sum 
rules fill]. Taking the 0(p*^) expressions for bi, 62 as input, he finds that at 
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0{p^) three of the bi are only weakly sensitive to the choice of the energy 
separating high- and low-energy components in the sum rules. The resulting 



values (for the interpretation of the errors we refer to [|TT|) 

63 = (-2.55 ±0.20) • 10-3 1 

64= (4.55 ± 0.15) ■ 10-3 > (5.11) 
be = (0.92 ± 0.03) ■ 10~^ ^ 

are consistent with our set II and with the values ( ^.10|) of Knecht et al. [|18 . 



5.3 Threshold parameters and phase shifts 

To compare the theoretical amplitudes with data on tttt scattering, one expands the 
combinations with definite isospin in the s-channel 

r°(s,t) = 3A{s,t,u) + A{t,u,s) + A{u,s,t) 

T\s,t) = A{t,u,s) - A{u,s,t) (5.12) 

T^is,t) = A{t,u,s) + A{u,s,t) 

into partial waves, 

00 

T\s,t) = 327r^(2/±l)P,(cos^)t/(s) . (5.13) 

1=0 

Unitarity implies that in the elastic regionQ 4 < s < 16 the partial wave amplitudes 
tj are described by real phase shifts 5/, 

tf(.) = (-^)'''l{e-^f(^)-l}. (5.14) 

The behaviour of the partial waves near threshold is of the form 

Retf(.)=g>f±g26f±0(g^)}, (5.15) 

with q the center-of-mass three-momentum of the pions, i.e. s = 4(1 ± /Ml). 
The threshold parameters a/^, b\ are referred to as scattering lengths and slope 
parameters, respectively. 

In table |, we compare the threshold parameters of the low partial waves with 
experimental results. We emphasize once more that we do not consider the values 
for either set I or set II as the definitive results of 0{p^). To arrive at such results, 



a more detailed analysis based on the Roy-equation approach is under way HT 
With this caveat in mind, we notice that the /i, I2 from the D-wave scattering 
lengths (set II) tend to improve the agreement with experiment with the possible 
exception of the S'-wave scattering lengths Oq, Oq. We comment on the experimental 
determinations of scattering lengths below. We also note that the ambiguity in scale 
or in the choice of the lagrangian of O(p^) affects aj] in the third digit only. 



^We recall our normalization of the Mandelstam variables in ( 4.13 ). 
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Table 4: Threshold parameters in units of Mt^+. The values of 0{p^) and 0{p^) (set 
I) correspond to h = -1.7, I2 = 6.1, [3 = 2.9, U = 4.3. For 0(/)_(set II), h = -1.5, 
I2 = 4.5 were extracted from the D-wave scattering lengths U unchanged). The 
fourth column shows the values obtained taking only the contributions from the fc, 
ll| defined in App. D. We set /i = 1 GeV and take the r[(l GeV) from The 
experimental values are from Ref. PU . 



0{p^) set I 



0(/) 
set II 



experiment 






26 ± 





05 





25 ± 





03 





28 ± 





12 





82 ± 





08 





38 ± 





02 



-10 al 
-10 bl 
10 a\ 
10%\ 
lO^aO 
lO^a^ 



0.16 
0.18 
0.45 
0.91 
0.30 









0.20 
0.25 
0.42 
0.73 
0.37 
0.48 
0.18 
0.21 



0.213 
0.279 
0.407 
0.69 
0.40 
0.78 
0.30 
0.34 



0.217 
0.275 
0.413 
0.72 
0.40 
0.79 
0.27 
0.23 



0.206 
0.249 
0.443 
0.80 
0.38 
0.54 
input 
input 



0.17 ±0.03 
0.13 ±0.30 



In Figs. H, P we plot the phase shift difference Sq — S} and the 1 = 2 S'-wave phase 
shift 6q as functions of the center-of-mass energy and compare with the available 
low-energy data. The phase shifts correspond to the amplitudes of O(p^), O(p^) 
and 0{p^) (for both sets I and II), respectively. The two-loop phase shifts describe 
the data quite well for both sets of U, with a small preference for set I. The 
1 = 2 S'-wave, on the other hand, seems to prefer set II. 

There is an interesting lesson we can draw from our amplitude concerning the 
extraction of threshold parameters from phase shift data. Looking at Fig. ^ we 
observe that the p^ amplitude (set I) describes the data almost perfectly. On the 
other hand, the value for the scattering length Og in table || is 0.217 for set I, quite a 
bit smaller than the 0.26 from Ref. 0] and even smaller than the 0.28 from Ref. |^ 



To put these differences into perspective, we consider the effective range approx- 
imation for S'-wave scattering 



gcot5o' = ^±kV, (5.16) 



with Tq the effective range. In terms of the threshold parameters defined in ( ^.15 ) 
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tree 

one loop - set I 
two loops - set I 
two loops - set II 



275 



300 



325 



350 



375 



E„^(MeV) 



Figure 8: Phase shift difference Sq — 6} at O(p^), 0{p'^) and O(p^) (set I and II). 
The data points are from Ref. p2| . 



the effective range is given by 



(5.17) 



by (g2a(, q%{) 



A similar formula can be used for 5{ when we replace 

In Ref. the following approximation was used to extract the scattering length 
Oq from the measured phase shifts: 



sin2(5°-5i) 




s 



(5.18) 



In addition, a relation between and h attributed to Basdevant et al. was 
employed for a one-parameter fit to the data. 

In Fig. |lO|, we compare the phase shift difference 5q - 5\ from the full two-loop 
calculation (set I) with the effective range formula ( |5.16| ) and with the approximation 
( [5.18[ ) used by Rosselet et al. For the effective range approximation, we have 
used the threshold parameters corresponding to the CHPT amplitude (set I in table 
H). In the curve based on ( |5.18| ) we have also used On = 0.217 and the above 



mentioned relation |^ between and h. 
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Figure 9: Phase shift 6q at O(p^), O(p^) and O(p^) (set I and II). The data points 
are from Ref. |E3 . 



The obvious conclusion is that both the effective range approximation and the 
approximation used in Ref. deviate from CHPT to 0{p^) already at compara- 
tively low energies. Fig. ^ also shows why the scattering length obtained in Ref. |^ 
is bigger than the CHPT value. It is hardly necessary to emphasize that the phase 
shift from CHPT is superior to both approximations on all accounts. 



6 Conclusions 

In this paper we have presented the calculation of elastic pion-pion scattering to 
sixth order in the low-energy expansion of QCD. The main part has been devoted 



to explaining the technical aspects of the results given in [|T0 

We first developed the loop expansion and the renormalization procedure for 
the generating functional of Green functions at the one- and two-loop level in the 
A^-component 0^ theory We concentrated on issues that are relevant for the corre- 
sponding CHPT calculation. In particular, we have discussed the dependence on 
the renormalization scheme and the impact of the EOM. Using the EOM, the coun- 
terterms in 0'^ theory can be written in different ways. We demonstrated explicitly 
that different choices are equivalent in the sense that the differences can always be 
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Figure 10: Comparison of th e full two-loop phas e shift difference 6q — 6l with the 
effective range formula ( ^.161) and with Eq. (|5.18|) . 



absorbed into the higher-order parameters. We discussed a method for calculating 
two-loop diagrams by treating the renormalized one-loop amplitude as a nonlocal 
vertex. This approach turns out to be especially useful for the more involved CHPT 
calculation. 

For the case of CHPT, we then performed the corresponding calculations for 
the pion mass, the pion decay constant and the tttt scattering amplitude. The 
renormalization of these quantities to 0{p^) was discussed in detail for both minimal 
subtraction and for the modified minimal subtraction scheme that we actually used. 
In analogy to 0'^ theory, the role of the EOM in relating different forms of the chiral 
lagrangian was discussed. We also commented on the Laurent expansion in ci — 4 of 
the coupling constants of O^p'^). 

We determined the complete dependence of M^r, F^r and the scattering amplitude 



on the quark masses. As in previous full two-loop calculations [|^, 0, this 
dependence is both of theoretical interest and of numerical relevance. The analytical 
expressions for the scattering amplitude and for the threshold parameters were given. 

For the numerical analysis of vrvr scattering, we derived estimates of the low- 
energy constants of 0{p^) on the basis of meson resonance exchange including kaon 
and eta contributions. These estimates turn out to be dominated by vector meson 
contributions. Although there are considerable uncertainties in this simple version 
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of resonance saturation, the overall size of the constants of 0{p^) is such that many 
quantities are relatively immune to those uncertainties. This is especially the case 
for the S-waves. However, for the higher partial waves and for a systematic error 
analysis in general a more sophisticated approach is needed. Such an analysis based 
on Roy equations is under way |4l|]. A similar approach has already been used in 
the dispersive treatment of the scattering amplitude to two-loop order [T 



The major uncertainties for a numerical analysis reside in the low-energy con- 
stants of 0{p'^). To illustrate this uncertainty, we have presented results both for 



the standard set of those constants 0, |2T|] and for a second set where /i, I2 are 
determined from the Z^-wave threshold parameters to 0{p^) as given in this paper. 
In this case, even the S'-waves are sensitive to which set of constants is used. This 
makes it all the more necessary to perform an analysis where all the constants are 
determined on the basis of 0{p^) quantities that include, in particular, quark mass 
effects. 

We have also presented numerical results for the phase shifts including the com- 
bination that can be measured directly in K^^ decays. For this particular example, 
we demonstrated that approximations like the effective range expansion may deviate 
from the CHPT predictions already at relatively low energies. We have emphasized 
the problems of extracting threshold parameters on the basis of such approximations 
from phase shift data. 

The results obtained here for the tttt scattering amplitude to 0{p^) will be used 
for a systematic analysis |^ of the available experimental data. In addition, they 



will be the basis for significant tests of QCD at low energies together with forth- 
coming measurements of vrvr scattering near threshold in K^a decays and in pionium 
decay P6|. 
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A Fish and sunset 

A.l Notation 

We use the notation 




together with 



diagrams 
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Furthermore, as only one mass parameter occurs in this appendix, we put everywhere 
= 1. It is straightforward to supplement the quantities below with the relevant 

mass factors. 

The loop function J(s) in d dimensions is 
/ 1 \ 



J{s) 



\{l-ll){l-{h-pY)l 
- 74^r(-^) £ dx (1 - sx{l - x)r ; p'-s. (A.l) 



(47r) 

For J one has the dispersive representation 

[dx 



J ^ f _\dx\_ _ ^ < w < , (A.2) 
Ja X — s 



where the measure is 



(47r)2+^r(3/2 + w) \4j ' 4 
It is often convenient to split off the divergent part through 



J(s) = J(0) + J(s) , 

;r(- 

[dx] 



^(0) - (4;^^(-^.) , 



J{s) = s 



oo 



1 



167r2 



4 X[X — S) 

' ; a = (l-4/s)^ . (A.4) 



cr In h 2 

a + 1 . 



The constant J(0) is related to Tm, used in the text, via Tm — —J{0) taking into 
account that M = 1 in this appendix. 

A.2 The fish diagram 

The scalar integral for the fish diagram has the form 

// 4 ^ \\ 



, i=l , 



with 



L>3 = 1-^2 , D^ = l~{l2 + h-pif , 

Q = Pl + P2 , = P2 = 1 > = s ■ 



42 



We evaluate this integral in two steps: First, we identify that part of the function 
V{s) that stays finite as ci — > 4 by subtracting the divergent subintegral and by 
removing the remaining overall divergence. In the second step, we determine the 
finite part by evaluating its absorptive part along the cut in the variable s and by 
then constructing an analytic function that has the same absorptive part and the 
same behaviour at s ^ 0, — oo. 



Subtractions 

Integration over I2 generates the loop function J{i),i = {pi — li)^. We subtract 
the divergence in this subdiagram by using the decomposition ( [A .41 ) and obtain 

V{s) = J{s)J{0) + Vi{s) , 

Vi is still divergent at ci = 4. Using the Feynman parametrization 
1 2x2dxidx2 



f 

Jo 



DiD2{x-i) Jo [(xiDi + {l-Xi)D2)x2 + {x-i){l-X2] 



^3 ' 



integrating first over li and then by parts in Xi shows that Vi can be made finite by 
subtracting its value at s = 0, 

Vi{s) = Vi{0) + Vi{s) ■ lim— = 0. 

Vi(0) can be evaluated by methods similar to the ones used below for the sunset 
integral. 

The finite part Vi{s) 

To evaluate the absorptive part of Vi{s), we invoke unitarity [ ^ : 
lmV{s) = J(0)ImJ(s) + ImFi(s) 

= J dfi{l)dfi{l')6'^''\Q-l-l')J{t) , 

where d/j, is the Lorentz invariant measure in d dimensions, 

J2W+3J I 

Ml) = 2(2^)2.^3^0 , ^° = V 1 + , . (A.6) 
We write again J{t) = J(0) + Jit) and use 

^ ' ' dn{l)dn{l')S^''\Q - I - I') = 

7rr(3/2) fsY f 



(47r)2+"'r(3/2 + w) V4y 4 
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This expr ession agrees with the absorptive part of J{s) in d dimensions, see Eqs. 
Q). Therefore, at d = 4, one has 



ImVAs) 



I dii{l)dii{l')5^^\Q - I - l')J{i) 



(A.7) 



where the measure in ( |A.7D is now the ordinary four-dimensional one, obtained from 
( [A .61 ) by putting w = 0. For J we insert its exph cit ex pression ( [A.4| ) at c? = 4. It is 
convenient to work out the phase space integral ( |A.7|) in the center-of-mass frame 
where 



Q'^ = (v^,0) , p^ = (v^/2,0,0,p) ,p=^s/A-l. 

After integration over I', one is left with an integral over the three-momentum I. 
Using the remaining one-dimensional delta function, all integrals except the one 
over the angle 6{l,pi) can be done easily. The integral left over is proportional to 



dz 



ai In ^— + 2 

(Ti + 1 



'1 



p{l - z) 



After the change of variables (cti — l)/(o"i + 1) = f , we find 



lmVi{s) 



4(47r) 



3(T + In 



1 - a 



1 - a 



a sa 



1 + cr 



a 



[l-A/s 



(A.8) 



It remains to construct a function with the proper cut structure and the correct 
behaviour at s ^ 0, oo. We find 



VAs) 



:i67r 



2\2 



IT 



with 



/ = crln :;— h iira. 



3sa^ 6 



l + a 



In the notation of Ref. 10 , this is 



K,(s) 



Of course, exactly the same method can be applied to integrals with a more com- 
plicated numerator. 



A. 3 The sunset diagram 

For the sunset diagram integrals of the type 



[i-m-m-ip-h-kf] 



(A.9) 
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have to be calculated. Here we again focus on the scalar integral H, the remaining 
ones can be done analogously. For H we also present the procedure for the infinite 
part. Using the (i-dimensional dispersion representation ( [A.2|) for J, we get 



poo pi 

H = I [ds] I dx 
Ja Jo 



[x + — x) — — x) — If]"^ I 

The dHi integrals can be performed easily and the result is 

poo pi 

H{p^) = / [ds'] / dxF2[z2] (A. 10) 

J4 Jo 



where 



T{-uj) 
(4vr) 

-Z2(P^) = X + s'{l — x) — p'^x{l — x). (A. 11) 



We subtract and add the two first terms of the Taylor series of -F2[-22] in p"^ around 
= 1 and obtain for the finite part at c? = 4 



H{p')-H{l)-ip'-l)H'{l)= J^^ds'^l--J^ dx}C2{s',x-p'), 
where we have introduced the kernel 

IC2{s',x;p') = ^lhn{F2h(p^)]-F2h(l)]-(p^-l)F^h(l)]} 

1 j^^Z2{jQ_^ip'^ -I)x{l-X)\ 



The integral / dx K,2 could further be done in closed form - the result amounts to a 
twice subtracted one-loop self-energy integral with two propagators (with masses 1 
and s'). In the text, we use 



1 A 

-L / , , L ^ I , , , 2\ 



Hip') = j^^2 1 ds' f-^l dx }C2{s',x-p') , (A.12) 

which is finite at = 1. We also need H{1) and H'{1), where the poles at = 
are contained. For the evaluation of the infinite parts a recursion relation can be set 
up by partial integration in x in (A. 10). This method gives 



V ^ (47r)4+2^ ^ ^ l2 4 8 ^ ^ ^ 

In Ref. the evaluation of the sunset integral is discussed in the general mass 
case. 
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B Off— shell four— point function in d dimensions 

We display the four-point function to one loop in d dimensions. More precisely, we 
use the sigma model parametrization (|3.3|) and define 



dxi dx2 dx3 e-'(P^"i+P2"2-P3"«-f4"*) (O|T0Xa^i)0^a;2)0'(x3)0"'(x4) |0) 



-T 



lm;ik ; 



J. 2 2 2 2\ 



witl 



Pl+P2=P3+P4 , 

S= ipi+ P2f , t = {pi- P'if , U = {pi- Pif . 

The wave function renormalization constant Z is the one appearing in the two-point 
function ([4.2|). In the standard isospin decomposition 

T^"^-'^'{s,t,u;plplplpl) = 5'''S''^Ais,t,u;plplplpl)+ cycl. , 
the scattering amplitude is obtained by putting all momenta on the mass shell, 



A{s, t, u) = A{s, t, u- M^, M'^, M^, 



F2 



+ 0{p'). 



We find 



F'^A{s,t,u;pl,pl,pl,pl) = 

{s - M^) {f^ + {2s + t + u- SM^) Ji(s)} 

+ { [pIpI + pIpI - t{pl +pI- t)] Ji (t) 

+ [A13A24 - t(Ai3 - A24 + t)] Ut) 



-2t 



pI+pI-u] J3{t) + {P3,P4,t) (P4,P3,m)} 



+8hpiP2 ■ P3P4 + [piP3 ■ P2P4 + PlPi ■ P2P3] 



+ 



s+pI-pI-^M^/2 



M 



(B.l) 



where 

J,{s) 

Ms) 

Remarks: 



-J(s) 
1 



Ms) 



, ■J2 



s{d-l) 



M 



s{d-l] 



M' --] J{s 



-T, 



M 



and Aik= pi - pi 



(B.2) 



®In Sect. 4, we express the Mandelstam variables in units of the physical pion mass squared. 
For simplicity of notation, we use in this appendix the standard definition. 
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i) For on-shell momenta, the result agrees at d = 4 with the result given in . 

ii) For off-shell pions, the amplitude is not finite at d = 4, in contrast to the off- 
shell amplitude considered in 0| , where the four-point function of pseudoscalar 
densities was considered. In that case, there are additional contributions pro- 
portional to the low-energy constants I3 and /4 that remove the remaining 
divergences. 

iii) It is not surprising that the above amplitude is not finite off-shell: The con- 
struction given in Ref. only guarantees that the Green functions built from 
quark currents are ultraviolet finite - Green functions of pion fields are un- 
physical objects, even if they occur at intermediate steps of a calculation, as 
in the present context. 

iv) Finally, we mention that these divergences in the off-shell amplitude do not 
generate nonlocal singularities in the two-loop calculation. 



C Scattering lengths and effective ranges 



From the explicit expression for the scattering amplitude in equation (|4.15|) , it is 
straightforward to evaluate the threshold parameters a[ and bj . Using the definition 
( ^.15|) , we find 



a 







7M2 



1 + 



X 



° 327rF2 r ■ 7 

70 4 5 2 157r2 



49 + 5 61 + 12 62 + 48 63 + 32 64 



192 



63 



-— + IO61 + 2462 + 9663 + 6464 + 

1/D 7 



~ 47rF2r + 4 



281 
~9~ 



+ 462 + 48(63+ h) 



77489 41357r2 10 - 592 - 6448- j j j 

+ — 61 + 62 + ^3 + 688 64 + 288 65 - 32 be 

81 72 3 9 9 



M2 



167rF2 



1 — X 



2 + 61 + 16 64 



+ x' 



262 
~9~ 



227r2 



+ 4 61 + 64 64 



8ttF3 



1 - 



X 



97 



+ - 



X 



10591 



1 

1457r2 
12 



- 2 62 + 48 6. 



11^ 64- 32 - 752 - - 
+ y &i - Y &2 + Y 63 + — &4 + 3266 
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+x 



7003 1697r2 



2160 2160 



— 61 ho 

12 9 



152- 
45 ^ 



364- - 
^64 + 3266 
45 



al 



where 



1 



49 



4807r3F4 L 360 

67 1277r2 



+ &3 + &4 



2160 



29^ 19^ 28. 188. „^ 

H &i H 62 H ^3 H ^4 + 8 fog 

432 60 90 45 45 



M2 



X 



2,3,4 



167r^6i_2,3,4 and 65,6 



(C.l) 



(C.2) 



The coefficients bi are displayed in appendix 0. The terms between the last square 
brackets in the expressions for aj and bj generate contributions of 0{p^) - these are 
beyond the accuracy we aim at here. In order to keep the formulae as simple as 
possible, we nevertheless retain them. In our numerical results these are removed. 



D The constants 6, 



The quantities bi in Eqs. ( [4.14|) and ( [4.15| ) stand for 



7 113 



+ X2 



167r2 



9^9^ ^ 9 ^ 108 



17 



1 35091 



216 167r2 1296 J 



48 



+ ^ [4A;i + 28k2 - 6h + l^h] + [32 ll + Ul^-b II] 11-8 + | 



2 1 2 

62 = -8r,+2i:--L -- 

^ ^ 3 167r2 9 



+ X2 



Il67r2 



-24 /[ 



166 „. 



/n 18 /o /'i 

92 



203, 317 

L + 

54 



1 1789 



3456 167r2 432 



-i [54A;i + 62A;2 + 15A;3 + 10^:4] - [32 + 4 - 5 Q II + r^j , 



11 17 

2n + -Z--L - — 

^2^2 167r2l2 

311 1 7063- 
+ 



+ X2 



167r2 



9 ^ 3 ^ 3 ^ 216 



69 1 2 1 67r2 864 
1 5 



+ 2 [4 1[ + Q II + - [38h + 30A;2 - 3^:4] 
6 



2^^ 6^ 167r2 36 



+ ^4l^ 



91^92 9 4^ 216 



17 



1 1655' 



3456 167r2 2592. 



+ 2/2'/^-^[A:i + 4fc2 + fc4]+r^} 



167r2 



31 145 625 , 7 1 66029 

/o H L H 

6 ^ 36 ^ 288 864 167r2 20736 



107 



— Ki 
16 96 



k2 + r 



5 ; 



where 



167r2 



_^r_^r 257^ _1 1 11375 

18 ^ 36 ^^ 864 ^ 432 167r2 20736 



—ki -k2 + rl 

48 96 ^ 



(D.l) 



1 , Ml 
m ■ 



167r2 yu2 

h = (4/[(/i)-7,L)L; 71 = ^ , 72 = ^ , 73 = -^ , 74 = 2. (D.2) 

We have denoted by /[ (r[) the renormahzed, quark mass independent couphngs 

from £4 (Cq), with 01 a—^ = ^— . In the text, we also use the parameters /j, 

defined by 



7i 



327r2 



M2' 



(D.3) 



The scale dependence of r\ is fixed by the requirement = 0. 
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